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i Abstract. A state selected at random from the Hilbert space of a many-body system 

is overwhelmingly likely to exhibit highly non-classical correlations. For these typical 
states, half of the environment must be measured by an observer to determine the 
state of a given subsystem. The objectivity of classical reality — the fact that multiple 
observers can agree on the state of a subsystem after measuring just a small fraction 
^1 of its environment — implies that the correlations found in nature between macroscopic 

systems and their environments are very exceptional. Building on previous studies 
of quantum Darwinism showing that highly redundant branching states are produced 
ubiquitously during pure decoherence, we examine conditions needed for the creation 
of branching states and study their demise through many-body interactions. We show 
that even constrained dynamics can suppress redundancy to the values typical of 
random states on relaxation timescales, and prove that these results hold exactly in 
(N( the thermodynamic limit. 
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Hilbert space is a big place, exponentially larger than the arena of classical physics. 
The Hilbert space of macroscopic systems is dominated by states that have no classical 
counterparts. Yet the world observed by macroscopic observers exhibits powerful 
regularities that make it amenable to classical interpretations on a broad range of scales. 
How do we explain this? 

The answer, of course, is that Hilbert space is not sampled uniformly; rather, the 
initial state and the Hamiltonian governing evolution are both very special. Quantum 
Darwinism [U [2] is a framework for describing and quantifying what distinguishes quasi- 
classical states awash in the enormous sea of Hilbert space. 

Typical macroscopic observers do not directly interact with a system. Instead, they 
sample a (small) part of its environment in order to infer its state, using the environment 
as an information channel [3j. Thus, when we measure the position of a chair by looking 
at it, our eyes do not directly interact with the chair. By opening our eyes, we merely 
allow them (and hence, our neurons) to become correlated with some of the photons 
scattered by chair (and hence, its position). 

Consider a system S with Hilbert space of dimension Ds decohered by a multi- 
partite environment £ = ®f =1 £%^ where each £i has dimension Dg. To understand the 
perception of classicality by macroscopic observers, it is of great interest to understand 
the quantum mutual information between S and some subset of the environment (a 
fragment) T = ® ieF £%^ where F C {1, . . . , N}: 



Above, Hj^ and Hst are the respective individual and joint von Neumann entropies. 
We denote the size of the fragment by = \F\ = fN, where / = IF/TV is the fraction 
of £ contained in T . The mutual information averaged over all J 7 of a given fractional 
size / is written as 



When the global state of S£ is pure, one can show [4J that this function is non-decreasing 
and anti-symmetric about its value at / = x /2. 

In the absence of preferred initial states or dynamics, the natural question is: what 
is the typical amount of mutual information between S and J 7 , and how does it depend 
on the fractional size / of the fragment? To be quantitative, we use the Haar measure on 
the space of pure states in the global Hilbert space % = S®£ of dimension D = DgDg . 
(This is the natural, unique unitarily invariant measure on this space.) Page's formula 
for the Haar-average entropy of a subsystem [5j |6j [7] can be used to calculate [4] the 
average of T over T~L. If we hold / fixed, we find that lim N ^ OQ (T(f)) n = if / < x /2. 
In other words, for a randomly selected pure state in the global Hilbert space, an 
observer typically cannot learn anything about a system without sampling at least half 
its environment. States that deviate (even by exponentially small amounts) from this 
property occupy an exponentially small volume in Hilbert space [8] as N oo. (This is a 
consequence of the mathematical phenomenon known as the "concentration of measure" 
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Figure 1. We investigate an environment of 16 spins Si coupled to a single system 
qubit S with Hamiltonian and initial state given by Eqs. Q and ([5|. A fragment T 
is a subset of the whole environment S . The couplings di and rriij were selected from 
normal distributions with zero mean and standard deviations = 0.1 and <j m = 0.001. 
Crucially, the interactions between S and the Si are much stronger than those within 
S. That is, <id ^> cr m . 

in high- dimensional spaces [9], which can be thought of as an abstract law of large 
numbers.) 

It's natural to define the redundancy R$ as the number of distinct fragments in the 
environment that supply, up to an information deficit 5, the classical information about 
the state of the system. More precisely, R$ = l//<s, where f$ is the smallest fragment 
such that I(fs) > (1 — 5)H™*, and i^ ax is the maximum entropy of S. The dependence 
on S is typically [10J only logarithmic. At any given time, the redundancy is the 
measure of objectivity; it counts the number of observers who could each independently 
determine the approximate state of the system by interacting with disjoint fragments of 
the environment. As described in the previous paragraph, typical states in T~L will have 
1(f) ~ for / < 1/2 and, by symmetry, 1(f) ~ 2i/if ax for / > 1/2 , so R$ ~ 2 for any 
S. That is, half the environment must be captured to learn anything about S. These 
states are essentially non-redundant, and make up the vast bulk of Hilbert space. 

1. Dynamics 

But of course, we know that observers can find out quite a bit about a system by 
interacting with much less than half of its environment. This is because decoherence 
is ubiquitous in nature [HJ [121 E2> E] and redundancy is produced universally by 
decoherence in the absence of coupling between different parts of the environment [10J. 
However, realistic environments can have significant interactions between parts, so it's 
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important to study these interactions and their effect on redundancy. To see how high- 
redundancy states form through decoherence and how they can relax to a typical non- 
redundant state, we consider a model of a single qubit S (D$ = 2) monitored by an 
environment 8 = &)f =1 £i of TV spins (Ds = 2) 

with Hamiltonian 

i j,k 

where di are the system-environment couplings and are the intra-environment 
couplings. We take the initial state to be 

|*o> = ^=[|t> + U>]®[| + l> + |-l>]® JV - (5) 

For clarity, we denote the states of S with arrows (| f), | \)) and the states of the Si with 
signs (| + 1),| — 1)). (There are several ways to relax this model for greater generality, 
but they are unnecessary for elucidating the key ideas. We discuss generalizations at 
the end of this article.) 
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Figure 2. (Caption next page.) 



The rise and fall of redundancy in decoherence and quantum Darwinism 



5 



Figure 2. (Previous page.) We study the spin universe described in figure [l] (a) The 
redundancy R$ is the number of fragments of £ that provide, up to a fractional deficit 
S = 0.1, complete classical information about the system. The exact redundancy is 
supplemented by an estimate based on the linearly interpolated value of X(f) to guide 
the eye. This can be compared to R§ ~ 2, the redundancy of nearly all states in 
the global Hilbert space. The vertical dashed lines mark five time slices, (b) The 
mutual information X versus fragment size ^J 7 , and the entropy H$ of the system, 
at five time slices corresponding to different qualitative behavior, (c) The complete 
mutual information X versus both fragment size and time t. The five time slices 
are marked with thick black lines. Low correlation (t = 0.5) for small times means 
the environment "knows" very little about the system. Each spin added to T reveals a 
bit more about <S, resulting in the linear dependence of X. Decoherence (t = 2) sets 
in near = (\^Nad)~ 1 = 2.5. By that time, the density matrix of S is approximately 
a mixture of the two pointer states | t) an d | 1} singled out by the interaction 
Hamiltonian. Mutual information is still nearly linear in and redundancy is of 
order unity. Mixing within the environment can be neglected because t <C cr" 1 = 1000. 
Quantum Darwinism (t = 10) is characterized by a mutual information plot that 
rises quickly to the classical plateau; the first few spins in a fragment give essentially 
all classical information, and additional spins just confirm what is already known. The 
remaining quantum information (above the plateau) is still present in the global state 
but it is effectively inaccessible, in that it can only be recovered by an unrealistic 
observer accurately measuring the joint state of almost all of £ . After t ~ cr^ 1 = 10, 
only order unity spins are needed to determine the state of S no matter how large N 
is, so Rs ~ N. In the absence of the couplings between environment fragments 
this situation would persist forever. (For some environments, such as photons, this is 
indeed the case.) Relaxation (t = 150) occurs near t ~ r m = (\^Ncr rn )~ 1 = 250. 
Mixing within the environment entangles any given fragment's information about the 
system with the rest of the environment, reducing the usefulness of measurements on 
that fragment. The mutual information plateau is destroyed, so redundancy plummets. 
Equilibrium (t = 500) is reached for t ~ a" 1 = 1000, when the actions associated 
with interaction between individual spin pairs in the environment reach order unity. 
The mutual information plot takes the non-redundant form characteristic of a random 
states in the combined Hilbert space of S£. An observer can learn nothing about the 
system unless he samples almost half the environment. 

We use a numerical simulation (N = 16) to illustrate the build up of redundancy 
from the initial product state, and the subsequent transition to a typical non-redundant 
state. (See figure [TJ) The couplings are selected from a normal distribution of zero 
mean and respective standard deviations and a m . Our key assumption to produce 
a high-redundancy state will be that S is coupled to the Si more strongly than the Si 
are coupled to each other (o^ 3> a m ). This is an excellent approximations for many 
environments (e.g. a photon bath p~5],[l6], where effectively a m = 0) but not all (e.g. a 
gas of frequently colliding molecules). This is the only condition that physically selects 
S as distinguished from the Si. For brevity, we'll call the timeframe t <C cr" 1 the pure 
decoherence regime and t a" 1 the (intra-environmental) mixing regime. (We have 
set h = 1. In this article, we refer to interactions between spins within the environment 
as "mixing".) 
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In addition to the two timescales a^ 1 and a" 1 set by the typical size of the 
interaction terms, we also are interested in the times = (a/A^) -1 an d = 
(V^^m) -1 which scale with the size of the environment. Roughly, t > a^ 1 and 
t > a" 1 are times for which the actions associated with couplings between individual 
spins (including the system qubit) are appreciable. The earlier periods t > and t > r m 
are the times for which the collective action of the TV environment spins (on the the 
system and the environment itself, respectively) is strong. 

Figure [2] shows the rise and fall of redundancy in the environment for our model, 
as well as the quantum mutual information between S and J 7 as a function of fragment 
size IF. The maximum entropy of S is one bit: i^ ax = In 2. The system is decohered, 
Hs ~ i?^ ax , when the environment becomes fully entangled with it, I(f — 1) ~ 2i7^ ax , 
and this holds after t ~ tj. However, the mutual information does not form a plateau 
indicative of redundancy until t ~ a^ 1 . The plateau at T « i7^ ax corresponds to 
approximately complete classical information about S available in most fragments T for 
B F not near or iV. But once enough time passes for the mixing to become significant, 
t ~ r m , this structure is destroyed and the plot takes the form characteristic of typical 
non-redundant states. 

To better illustrate what is going on, the average entropy Hj?(f) is plotted in 
figure [3^. During pure decoherence, Hjr saturates at i^ ax for away from and 
N. However, once the mixing during relaxation becomes substantial, Hjr approaches its 
maximum values consistent with the dimensionality of T and the symmetry [Hjr = H s y 




Figure 3. (Caption next page.) 
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Figure 3. (Previous page.)We study the spin universe described in figure [l] (a) The 
average entropy Hjr of the fragment state pjr versus is the essential component of 
the mutual information, X = H$ + Hjr — Hst, for understanding the rise and fall of 
redundancy and how it relates to branches in the global state of SS. Five time slices 
are are marked with thick lines, (b) The three components of the mutual information 
(H s , and H S t) for each time slice. Note that Hjr^J 7 ) = H S t(N - IF) by the 
Schmidt decomposition and that — for all times after S is initially decohered — H$ is 
essentially equal to i^ ax = In 2. (c) The eigenvalues of the state pjr (which determine 
Hjr) for the same five time slices. The n-th largest value plotted is the average of the 
n-th eigenvalue of each choice of T . Low correlation (t = 0.5) exists before there are 
significant interactions, and there is just a single dominant branch corresponding to the 
initial product state. Decoherence (t = 2) produces two branches in the global state, 
one for each of the pointer states | t) and | \) of <S. The system has been decohered by 
the environment at this point, but only very large fragments are fully correlated with 
S. Observers measuring less than half of the environment will not be able to deduce the 
state of the system; as yet, there is no objectivity. The widely separated eigenvalues for 
$ T < N/2 imply that the global branch structure is not accessible to local observers. 
Quantum Darwinism (t = 10) is characterized by the fact that even small fragments 
T reveal the state of S — and hence which branch the observer is on. Since T only 
interacts with <S, there can only be two branches and the entropy Hjr is bounded by 
^max _ i n 2 By symmetry, the same is true for Hst- The tiny eigenvalues rising from 
below are the early indications of mixing. Relaxation (t = 150) causes the number 
of significant eigenvalues to expand because T now entangles with its complement 
F — ®i^F & m addition to S. The two branches (corresponding to the two dominant 
eigenvalues of pjr) are beginning to divide, so that knowing a small fragment of the 
original branch no longer suffices to specify the pointer state of its root. The entropies 
Hjr and Hst quickly exceed i^^ ax = In 2. Equilibrium (t = 500) follows. The state 
pjr approaches the maximally mixed matrix for < N/2, so the eigenvalues of pjr 
are clustered around 1/ dim (J 7 ) = 2~ T . The entropy Hjr approximately saturates its 
maximum, The global branch structure is destroyed and the composite system SS 
cannot be given a classical description. 



where T = ® i( £ F £i) imposed by the Schmidt decomposition: 

^max = Min [ttjr N + i _ Sjr] ln 2 , (6) 

H™* = Min[lF + 1, N - IF] In 2. (7) 

In figure |3]3, the eigenvalues for the corresponding state pjr are likewise plotted in 
both regimes. This shows the formation and destruction of branches characteristic 
of quantum Darwinism [TJ El [3j H] , and is suggestive of Everett's relative states pTTl IT8] . 
For pure decoherence, there are two dominant eigenvalues, corresponding to the entropy 
Hjr capped at i7if ax . As the mixing becomes important, the number of significant 
eigenvalues of pj=- quickly rises and pushes the entropy to its maximum. 
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Figure 4. The decoherence factor T, (10), for coupling constants di chosen from a 
centered normal distribution with standard deviation <Jd = 0.1. T can be expressed as a 
product of cosines that are statistically independent for times t ^> cr^ 1 . The statistical 
behavior of T is seen in this long-time plot, and it can be shown that fluctuations away 
from zero are exponentially suppressed in the thermodynamic limit. The dashed lines 

(17}. (a) TV = 8. (b) N = 12. The insets shows a 
16 would be too small to plot. 
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2. Branching 

We can develop a good intuition for this behavior by considering branches in the global 
state [HJ US] of SE. Suppose that at a given moment the state can be decomposed as 

m = Et*i^> = E (8) 

q=l q=l 

for some small number Q of orthogonal product state branches \ip q ). For t = 0, we 
can have Q — 1 since the initial state is a product state. In the decoherence regime 
(with approximate equality) we can have Q — 2, i.e. a generalized GHZ state [20J. 
But once the environment begins to mix, Q > 1. This gives a way for understanding 
the proliferation of eigenvalues in Hjr. For any choice of fragment J 7 , its entropy Hjr 
is bounded from above both by i7j£ ax ^ and by the entropy of the branch weights 
|7 g | 2 , because the Schmidt decomposition associated with the cut T-ST cannot have 
any more than Q branches. (See figure |3j) More precisely, the spectrum of the fragment 
state pjr cannot be more mixed than the probability distribution |7 g | 2 according to the 
majorization partial order [211 E2] for any choice of T . 

With this intuition in hand, we now derive the behavior seen in our model in the 
next two sections for large TV; mathematical details can be found in the Appendix. 

3. Pure decoherence 

In the pure decoherence regime, t <C cr" 1 , both decoherence [23j [24] and quantum 
Darwinism [3l [25j [26j [27] are well understood (even with Ds > 2). The single 
decoherence factor of the two-state system quantifies the suppression of the off-diagonal 
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terms of the density matrix pf with time: 

At) - \ ( ± r ? ) , (9) 




N 



r = Y[cos(2dit). (10) 

i=i 

The entropy of the two-dimensional state ps (|9j) is then 

H s = In 2 - r arctanh T - In Vl - T 2 (11) 

~ln2-^r 2 , (12) 

where the approximation is valid for small Y. The average mutual information between 
S and T is 



1(f) = H S + fl>(/) - ^(/) (13) 

«in2-^(r 2 + r3 ? -r^), (14) 

where r> = flieF cos(2dji), and F-= = Yl^ F cos(2e^). 

The short and long time limits are illuminating. For £ <C cr^ 1 and large iV, 
cos(2dit) « 1 - 2d 2 t 2 so 

r » exp(-2t 2 ^ d 2 ) - exp(-2tVJiV). (15) 

Therefore, the system is essentially decohered when t ~ r^, and i?^ ~ 1. The ensuing 
period < £ < a^ -1 exhibits quantum Darwinism. The system remains decohered but 
each spin in the environment continuously collects more and more information about 
the system. Consequently, the redundancy steadily rises because the number of spins 
that must be measured by an observer to determine the state of the system falls. This 
continues until t ~ cr^ 1 , when the phases associated with the action of the E{ on S are 
of order unity. At this point, the classical plateau of the mutual information congeals 
and Rs ~ N. 

We can be precise by looking at t 3> cr^ 1 , when the values of the cosines on the 
rhs of (10) will act as independent random variables [28j. The statistical behavior is 
described by the time-averaged expectation values 

«r» = o, (16) 
4 = «|rf » = 2-", (17) 

since ((cos 2 (2^i))) = 1/2. 

In other words, the decoherence factor Y has a Gaussian fall from unity for short 
times, and fluctuates around zero thereafter. This is illustrated in figure [4j The 
fluctuation of Y away from zero are exponentially suppressed, so fluctuations of Hs 
away from ifg 18 * = In 2 are similarly tiny. IV and Yy have the same behavior [with the 
respective replacements TV —> fN and TV — >■ (1 — f)N] so 

1(f) « In 2 - \ [2~ N + 2~f N - 2~^ N ] (18) 
->ln2 (19) 
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for / 7^ 0, 1 in the thermodynamic limit. This forms the robust classical plateau at 
1 = In 2 = H^. 

Although we concentrate here on the large time limit t 3> cr^ 1 for the sake of rigor, 
note that the plateau starts forming at t ~ and finishes at t ~ cr^ 1 . Indeed, even weak 
interactions lead to reliable redundancy [10J, a result that holds for higher dimensional 
subsystems. In particular, the ubiquitous real-life case of collisional decoherence through 
scattered light [T5l [IS] demonstrates how many weak correlations add up to huge 
redundancies. 



4. Mixing within the environment 

In the mixing regime, t a" 1 , interactions within the environment force distinct records 
about S stored in the Si to intermingle, making it more difficult on average to determine 
the state of S by sampling a given fragment T . For large times, the mutual information 
between S and a typical T is nearly zero unless / > x /2, i.e. an observer is unable 
to tell anything at all about the system until he makes a measurement on almost half 
the environment. Although the same amount of entanglement and information exists 
between S and S regardless of mixing within the environment, the mixing spreads 
this information globally, rendering it locally inaccessible. Information about S is no 
longer confined to the subsystems of £, but is stored in the correlations between them. 
Similarly, one learns nothing about whether or not a pair of playing cards are the same 
suit by looking at just one card. 

To see this analytically, we now show that pjr will tend to the maximally mixed 
state p^ = I/2^ N for large times. First, note that pj=- agrees with p^ on the diagonal in 
the z-basis |f), where \r) = (£)j eF \rj), rj = ±1, is a state of T specified by the vector 
f. The off-diagonal elements of pj=- are suppressed by the factors 



cos I tV(m jfe + m kj )(rj - r') 



jeF 



(20) 



k£F 

which are analogous to T. For t 3> a" 1 , the cosines will act like independent random 
variables and tend to cancel. To be specific, ((|A(^_p)| 2 )) = 2 - ( 1- -^ iV for r ^ r 1 . 

For large times, one can show that the chance of an exponentially small fluctuation 
in pjr away from the maximally mixed state becomes exponentially unlikely in the 
thermodynamic limit: 

Pr[r(p^p-)>e-^]<e- 2 ^, (21) 

where T = T{p J r ) p^) = \\pjr — p™\\i/2 is the trace distance and k is a strictly positive 
constant for / < 1/2. It is in this sense that p? approaches the maximally mixed state 
for / < 1/2. The Fannes-Audenaert inequality [29, 30j then implies that exponentially 
tiny fluctuations in Hjr are likewise exponentially unlikely over large times. In that 
sense we say that 

Hr^H\p?] = fNhi2 (22) 
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as TV oc for all T with / < 1/2. With only a minor modification, the same argument 



can be applied to pst to show Hst —± (fN+1) In 2. We know from ( 12 ) that H$ In 2, 
so 

Is-.? = H s + H T -H S T^0 (23) 

for all fragments satisfying / < 1 /2. Since Hg = H$ —> In 2, we get X(f = 1) —> 2 In 2, 
and so by the anti-symmetry we know 1(f) 2 In 2 for / > 1/2. 

This explains the persistent step-function shape of the mutual information for large 
times as plotted in figure [2j This is the same form of the mutual information obtained 
with overwhelming probability by a state selected randomly from the global Hilbert 
space. 



5. Discussion 



Decoherence [T4l fT3l [T2] is crucial for understanding how classicality can arise in a 
purely quantum universe. However, concentrating on individual systems (even while 
accounting for their interaction with the environment) leaves much to understand about 
global states. 

Quantum Darwinism has sharpened the vague idea that, based on our everyday 
observation of the effectiveness of our classical view of the world, there must be 
something very special about quasi-classical global states. Hilbert space is dominated by 
non-redundant states, and these are not consistent with the high redundancy observers 
take for granted when they extrapolate an independent reality based on local interactions 
with the immediate environment. 

Quantum Darwinism shows how high redundancy can arise from decoherence. 
However, in many-body systems branching states with large redundancy cannot last 
forever. The average mutual information 1(f) approximates a step-function for almost 
all the states in Hilbert space, so sampling ergodically produces such states with 
near certainty. Therefore, relaxing to equilibrium necessarily means the destruction 
of redundancy. 

If desired, our model can be generalized. "Unbalanced" initial states of the system 
[T6] or of the environment [27], such as (2| f) +i\ \))/V5 : do not change the qualitative 
results. The mutual information plateau will form lower at i^ ax < In 2 to agree with 
the maximum entropy of the system, and the limiting state p^ will change, but the 
factors A(f?_f*) controlling fluctuations in pj=- away from p^ will still be exponentially 
suppressed. The general unbalanced case is handled in the Appendix. 

We emphasize that the commuting nature of the interactions is very natural; the 
interaction terms between macroscopic objects (scattering) are almost always diagonal 
in position, a fact that can be traced back to real-world Hamiltonians. Adding a self- 
Hamiltonian for S or the Si diagonal in the z-basis will not change any of our information 
theoretic results, since all the relevant density matrix spectra will be the same. Self- 
Hamiltonians for S that do not commute with (PD) partially inhibit decoherence itself 
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[3TL [T2l [32] . but will not stop the information mixing in the environment. In general, 
system self-Hamiltonians that do not commute with the system-environment interaction 
are necessary to produce the repeated branching that occurs in nature. For example, 
the rate of diffusion for the quantum random walk of an object decohered by collisions 
with a gas is set by the size of the self-Hamiltonian p 2 /2m relative to the strength of the 
scattering [33]. An enticing subject for future research would be the analysis of quantum 
Darwinism in the case of repeated branchings due to a non-commuting self-Hamiltonian, 
and the dependency of redundancy on the rate of branching. In particular, we expect 
strong connections [SUES, 36j with the quantum trajectories [37] and consistent histories 
[381 139] formalisms. 

Our simple model has highlighted how important the relative strengths of couplings 
are for the distinction between system and environment, and the development of 
redundancy. Indeed, coupling strength is the only thing here that distinguished 
the system from the environment. If we had not assumed that the mixing within 
the environment was slower than the decoherence of the system, there would be no 
intermediate timespan a^ 1 <C t <C cr" 1 and the mixing would destroy redundancy 
before it had a chance to develop. 

Such mixing would seem unimportant when studying the decoherence of a system of 
a-priori importance, but it's illuminating for understanding what distinguishes certain 
degrees of freedom in nature as preferred. A large molecule localized through collisional 
decoherence by photons is immersed in an environment with insignificant mixing [40J, 
and so is recorded redundantly p~5j HE] , but a lone argon atom in a dense nitrogen gas is 
not. Whether an essentially unique quasi-classical realm jHJHS] can be identified from 
such principles is a deep, open question j43j[44] about the quantum-classical transition. 
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Appendix 

Here we discuss decoherence factors in the thermodynamic (N 00) and large time 
(t 00) limit. Recall that our model consists of a single qubit S monitored by an 
environment S = &)f =1 Si of N spins 

n = S®E^C 2 ® {C 2 ) m (A.l) 

with Hamiltonian 

i j,k 
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where d{ are the system- environment couplings and and the environment- 

environment couplings. The initial state is 

l*o> = t) + h\ ;>] ® Kil + 1) + - 1>]®" (A.3) 

where |/3^| 2 + |/3^| 2 = |<^+i| 2 + |<^_i| 2 = 1. Let T = ® ieF £i be a fragment of £, where 
F C {l,...,Aff and ^J 7 = |F| = /AT, < / < 1. The complement fragment is 

Let us break up the evolution into commuting unitaries 

U = e -H* = u 5J .u 5 _u^u^_u_ 5 (A . 4) 

labeled by the subsystems they couple, e.g. 



Up? = exp 



The single decoherence factor of the two-state system quantifies the suppression of the 
off-diagonal terms of the density matrix p s with time: 



where 



N 



T(t) = H [\a +1 \ 2 e~ mdi + \a^\ 2 e 2itdi ] . (A.7) 

i=l 

We are interested in the statistical behavior of this term for times large compared to 
the d i) especially for large values of N . For any function /i(i), we can define a random 
variable Z over a rigorous probability space through the cumulative distribution function 

F z (z) = P[Z >z}= lim h{t e [0,TMt) > z}, (A.8) 

T-»oo 1 

provided the limit exists. (Here, A is the Lebesgue measure). To be suggestive, we 
can denote expectation values over long times constructed with such a random variable 
using the time-dependent function: ((/x 2 )), ((log(/x))), etc. A result of the theory of almost 
periodic functions [451 146] is that random variables defined in this way from periodic 
functions of time are statistically independent if their periods are linearly independent 
over the rationals [28J. Unless the di are chosen to be exactly linearly dependent, this 
means that 



«r» =0, (A.9) 

N 

«|r| 2 » = n<<l ae " 2 ^ + (! - o)e md ^\ 2 » = [a 2 + (1 - a) 2 ] N , (A.10) 
i=i 

since ((e~ 4ltdi )) = 0. We have defined the probability a = \a+i\ 2 and note that 
y 2 <a 2 + (l-a) 2 < 1. 

Thus, so long as the environment isn't an exact eigenstate of the Hamiltonian 
(a 7^ 0,1), fluctuations of the decoherence factor T away from zero (as measured by 
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the variance) are exponentially suppressed in the thermodynamic limit. (For physical 
intuition about these results, see [24] .) This sends 

(A.ll) 



where H 2 [x] = - 
Of course, #™ ax 



xlnx — (1 — x) ln(l — x) is the binary entropy function and b = |/3^| 2 . 
< In 2, with equality iff b = \. 
We can quickly extend this to a statement about quantum Darwinism in the case 
of pure decoherence (rrijk negligible): 

*V (Po ® Po ) 
^t)(^t| + |/3;| 2 |^l)(^l| 



Tr 



ST 

2| 



where \T t) = ®j^ F { a ^ e 



-idnt\ 



(A.12) 
(A.13) 

1) + a_ie +ldjt \ — 1)) is the pure state of T conditional 
on the system being up, and likewise for \T \). The decoherence factor of the rank-2 
matrix pjr is Yjr = (J 7 I \ T t) and 

iv = n i a e ~ mdj + - a y 2itdj ] , ( a - 14 ) 

=0, (A.15) 
((\r^))=[a* + (l-ay] fN . (A.16) 

For fixed / > 0, fluctuation in IV will be exponentially suppressed in N, just like V. 
With small I>, 

2 



I 



it max 
^5 



where 



X(b) 



^1 - 46(1 - 6)(1 

|p 12 

46(1 - 6) arctanh (1-26) 



ir 



1) 



(A.17) 
(A.18) 

(A.19) 



1 - 26 

and < x(b) < 1 with x(b) = 1 iff 6 = Y2. Therefore, _£/jf quickly approaches H$ 
Further, because the global state is pure, we know Hsf = Hy and so 



l-s-.T = H s + B.J — H S jr 



rrmax 



X(b) 



iri 2 + ir 



(A.20) 
(A.21) 



For / away from and 1, all three decoherence factors are exponentially suppressed in 
N . This is the origin of the robust plateau on the plot of average mutual information. 
One can show [16j this means the redundancy grows linearly with N. 

Now we will extend this result to determine the statistical behavior of Hjr and 
Hsf when the interactions within the environment are not negligible. This will let 
us show that for times large compared the couplings the states of T and ST 
become maximally mixed subject to constraints of the initial conditions. First, under 
the evolution of H , the state pj=- of the fragment is unitarily equivalent to 



(A.22) 



The rise and fall of redundancy in decoherence and quantum Darwinism 



15 



A bit of algebra gives 

(r\pA^) = [fee-^^F^-^-rj) + ^ _ b ^ e +itJ2 jeF d j (r j -r>)'] JJ a r . a * 

where \r) = (££L GF | r ^ ) , = ±1, is a state of T specified by the vector r. Above. 



(A.23) 



JJ aexp -it^2(m jk + m kj )(rj - r'j) 



k(£F 



(A.24) 
(A.25) 



+ (1 - a) exp I ^2(m jk + m kj )(rj - r£ 

We want to show that the entropy Hjr of pjr approaches its maximum value fNH 2 [a] 
for / < !/2 by bounding the difference between p? and the limiting state pj?\ 

= ® [a\ + 1)(+1| + (1 " a) | - 1)(-1|] • (A.26) 



ieF 



First, we will assume the case of a balanced initial environmental state, a = 1 /2. The 
Hilbert-Schmidt norm of the difference is 



i^-^iils = EEi^i r ")i s 



r r'^r 



< 



^EEiv- 



4fN 



■r") | 



(A.27) 
(A.28) 



r r'^r 

Now, we want to bound fluctuations of Hjr away from its limiting value H™, the entropy 
of pj?. To do this, we use Audenaert's optimal refinement [30j of Fannes' inequality [29] 
governing the continuity of the von Neumann entropy. For any two density matrices p\ 
and p 2 with trace norm distance T = T(pi,p 2 ) = \\Pi ~ P2II1/2, the difference in their 
entropies AH is bounded as 

\AH\ < T\n(D - 1) + H 2 [T] (A.29) 

where D is the dimension of the matrices. We will also use the bound between the trace 
norm distance and the Hilbert-Schmidt norm for Hermitian matrices, | \p\ |i < VD\ \p\ |hs, 
to get 



T < l -V2^\\~ Pr - ~p™\ 



Ihs- (A.30) 

Now we consider the likelihood of fluctuations in T bigger than an arbitrary T : 

2T 



P[T >T }<P 



|hs > 



2/JV/2 



\Pr- PtWIs > 



4T 2 



2f N 



By the definition of an expectation value, we know 



5.F - P^Hls > 



4T[ 
2f N 



ATr 



2f N 



<«IIp>-p?IIhs) 



(A.31) 



(A.32) 
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And, from (A. 27), we know 



i^-/^ii!s>><^E£«i a 



AfN 0(1-/ 



< 



4/w Z-^^2( 1 -/) 7V 

l 



(A.33) 
(A.34) 
(A.35) 



We calculated 



\^(r-r') 



2(i-/)JV • 

| 2 )) exactly as we did ((r)). If r 7^ f* ', the the sums over 



j in (A. 24) are non-empty and — assuming the rrijk aren't specially chosen to be 
linearly dependent over the rationals — each fc-indexed term in the product of (A. 24) 
is statistically independent. 

Combining flA3I| ), §AJS2§ , and (|A35[, we find 

2 (2/-l)iV 



P[T > To] < 



4T 2 



(A.36) 



So if / < V2, we can choose T = 2"( 1 /2-/W2 so that both T and P[T > T ] are 
suppressed: 

P[T > 2"( 1 /2-/)iV/2 ] < 2 -(l/2-/)iV_ (A _ 37) 

In other words, as we take the size of the environment N to infinity, exponentially tiny 
fluctuations in the trace norm distance T = T(pp, p™) become exponentially unlikely. 
It is in this sense that we say 

Pt ~Pt (A.38) 

up to unitary equivalence. We can slightly relax the Fannes-Audenaert inequality to 
make it a little more transparent: 



|A#| <T ln(D-l) + # 2 [T ] 



l + /AHn2 + ln 



1 



(A.39) 
(A.40) 



So likewise for the entropy Hjr, exponentially tiny fluctuations are exponentially unlikely 
for large N. It is in this sense that we say 

H T ^H? = fNhi2 (A.41) 

for / < 1/2. 

With only a minor modification, the same argument can be applied to pst to show 



tt v rj cxo Zjoo 



it max 



(A.42) 



We know from Q that # 5 -> #s a3 S so J 5: ^ = H s + H r - H ST ->• for / < 1/2. 
Since = ; we get X^.jr — >■ 2i/™ ax for / = 1, so by the anti-symmetry we know 
2s : jf — >■ 2i7™ ax for / > 1/2. This gives exactly the step-function-shaped curve of a typical 
non-redundant state, so that R§ ~ 2, independent of S. 
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a 



Figure Al. Our argument is valid for / < /*, where /* is a function of how balanced 
the initials state of the environment is, as parameterized by a = |«+i| 2 . Only a < 1 /2 
is shown because — a) = /*(a). The colors denote the value of /* without (black) 
and with (red) Schumacher compression. 



If we directly extend this proof to the unbalanced case, a ^ 1 /2 ) we make the 
replacement 

1 r \2lO--f)N 



[a + (l-a) 2 ]' 



(A.43) 



in (A.35), but the factor of 2^ N from (A.30) is unchanged. This means we can show 
Hjr -> ifj? = fNH 2 [a] 



only for / < /* where 



/* 



(A.44) 



(A.45) 



log 2 [a 2 + (l-a) 2 ]J • 
Note that < /* < x /2, and /* = x /2 iff a — x /2. We can use Schumacher compression 
[13 HE] to improve /*. Take the typical sequence [19] of eigenvalues Af = YijeF \ a rj\ 2 °f 

Ss = {f I e-fWiM+S) < Xr < e -/W"M)}, (A.46) 

and define the typical subspace As C J 7 , S > 0, as the subspace corresponding to those 
eigenvalues 

A s = span{|A f ) | r £ S^} (A. 47) 

with projector II^ = J2? eS§ Then define P$ = = TiUsp^Us and the 

(normalized) density matrices 

fa = ^rUsprlls, (A.48) 



1 



■n^n*. 



Use the triangle inequality to bound 

ll^-^lli < ii^-^iii + ii^-^iii + ii^-P~iii- 



(A.49) 
(A.50) 
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The norms \\pjr — fjjr\\i and \\fjj? — pj?\\i can be handled with the close relationship 
between the fidelity and the trace distance: 1 — F(pi,p 2 ) < D(pi,p 2 ) < y/l — p 2 ) 2 - 
Using Hoeffding's inequality [50] we can show that these norms are suppressed 
exponentially in N. Now we just bound \\fjjr — ff^\\\ using the Hilber-Schmidt norm 
where, importantly, fjj- and fj^ live in a subspace of dimension e^ NH2 ^ rather than 2^. 
This gives an improved range for the applicability of our argument: / < /* where 

H 2 [a] 



u 



1 



(A.51) 



log 2 [a 2 + (l-a) 2 ]J • 

The improvement on /* is depicted in figure Al. This turns out to be the best we can do 
using the bound (A. 16). It's possible to construct a counter-example matrix pjr when 
/ > /* that satisfies ( |A.lo| but has limiting entropy Hf = (l-/)Aln[l/(a 2 +(l-a) 2 )] < 
Hf. 

So, in the case that a ^ 1 /2, we are only able to prove that Hjr —> Hj? = fN In 2 and 

1(f) for / < This means the redundancy can be bounded only by R$ < Iff*- 

Now, f* is of order unity unless the initial state of the environmental spins are nearly 

eigenstates of the interaction Hamiltonian, so this is still a very strong upper bound on 

the redundancy. In contract, R$ grows linearly with TV for a branching state. 
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